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FOREWORD 


The  study  reported  herein  was  conducted  at  the  U.  S.  Army  Engineer 
Waterv/ays  Experiment  Station  (WES)  during  the  period  1970-1971.  Guid¬ 
ance  and  funds  for  this  work  were  provided  by  the  Office,  Chief  of  Engi¬ 
neers,  under  Military  Engineering  Research  Project  No.  4A061102B52E-01, 
and  by  the  U.  S.  Army  Engineer  Topographic  Laboratories  under  Military 
Geographic  Analysis  Project  No.  4A062 1 12 A8 54 -02 . 

The  tests  were  performed  by  personnel  of  the  Terrain  Analysis 
Branch  (TAB),  Mobility  and  Environmental  Systems  Laboratory,  under  the 
general  supervision  of  Messrs.  W.  G.  Shockley  and  S.  J.  Knight,  Chief 
and  Assistant  Chief,  respectively,  of  the  Mobility  and  Environmental 
Systems  Laboratory;  W.  E.  Grabau,  Chief,  TAB;  and  A.  N.  Williamson, 
Project  Manager.  Mr.  J.  R.  Lundien  directed  the  study  with  the  assis¬ 
tance  of  Dr.  H.  Nikodem.  Persons  actively  engaged  in  the  program  in¬ 
cluded  Messrs.  E.  A.  Baylot,  J.  R.  Lundien,  P.  A.  Smith,  and  J.  H. 

Meeks  and  Miss  aret  H.  Smith  of  TAB,  and  Mr.  G.  H.  Jetton  of  the 
Instrumentation  Ser.lces  Division,  WES.  This  report  was  written  by 
Mr.  Lundien. 

Special  acknowledgment  is  made  to  representatives  of  the  Missis¬ 
sippi  State  Highway  Department  for  their  coorporation  in  providing  test 
sites  and  highway  design  information  for  this  test  program. 

COL  Ernest  D.  Peixotto,  CE,  was  Director  of  the  WES  during  the 
conduct  of  this  study  and  preparation  of  this  report.  Mr.  F.  R. 

Brown  was  Technical  Director. 
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NOTATION 


Velocity  of  electromagnetic  wave  through  free  space,  m'sec 
Layer  thickness,  m 
Base  of  natural  logarithm 
Electric  field  intensity,  v/m 
Incident  electric  field  intensity,  v/m 
Reference  electric  field  intensity,  v/m 
Reflected  electric  field  intensity,  v/m 
Freauency,  Hz 
\£ 

Harmonic  number 

Reflected  power  from  the  reference  standard,  w 
Reflected  power  from  the  test  site,  w 

Transmitted  power  with  the  reference  standard  in  place,  w 
Transmitted  power  with  the  test  site  exposed,  w 
Electric  wave  reflectance 
Power  reflection  coefficient 

Power  reflection  coefficient  for  perpendicular  polarization 
Power  reflection  coefficient  for  parallel  polarization 
Time ,  sec 

Oscillation  period,  Hz 
Distance,  m 
Impedance,  ohms 

Attenuation  constant,  m” 

-1 

Phase  constant,  m 
Propagation  factor,  m'^ 

Dielectric  constant,  F/rn 
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sq  Dielectric  constant  of  free  space,  F/m 

e  Relative  dielectric  constant,  dimensionless 

r 

0  Angle  of  incidence,  deg 
X  Wavelength,  m 
ix  Magnetic  permeability,  H/m 

|xq  Magnetic  permeability  of  free  space,  H/m 

Relative  magnetic  permeability,  dimensionless 
Conductivity,  mhos/m 
co  Circular  frequency,  radians /sec 
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CONVERSION  FACTORS,  BRITISH  TO  METRIC  AND 
METRIC  TO  BRITISH  UNITS  OF  MEASUREMENT 
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British  units  of  measurement  used  in  this  report  can  he  converted  to 
metric  units  as  follows: 


_ Multiply 

Fahrenheit  degrees 
inches 


By  _ To  Obtain _ 

5/9  Celsius  or  Kelvin  degrees* 

2.54  centimeters 


Metric  units  can  be  converted  to  British  units  as  follows: 


_ Multiply _ 

meters 

grams  per  cubic  centimeter 

centimeters 

millimeters 


By  _ To  Obtain 

3.2808  feet 

62.43  pounds  per  cubic  foot 

0.3937  inches 

0.03937  inches 


*  To  obtain  Celsius  (c)  temperature  readings  from  Fahrenheit  (F)  read¬ 
ings,  use  the  following  formula:  C  =  (5/9)  (F  -  32).  To  obtain 
Kelvin  (K)  readings,  use:  K  =  (5/9) (F  -  32)  +  273.15. 
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DETERMINING  PRESENCE,  THICKNESS,  AMD  ELECTRICAL  PROPERTIES 
OF  STRATIFIED  MEDIA  USING  SWEPT-FREQUENCY  RADAR 

PART  I:  INTRODUCTION 

✓ 

Background 

1.  Quantitative  information  on  surface  and  subsurface  materials 

is  used  for  a  variety  of  purposes:  vehicle  performance  or  foundation  de¬ 
sign  strength  can  be  predicted  if  the  properties  of  the  soil  or  rock  in¬ 
volved  are  known;  and  the  carrying  capacities  of  pavements  and  other 
artificial  surfaces  can  be  evaluated  if  the  strength  of  the  medium,  layer 
thicknesses,  foundation  properties,  etc.,  can  be  determined.  The  nature 
.and  extent  of  the  sampling  and  testing  program  necessary  to  acquire  this 
.  information  must  be  determined  by  means  of  the  following  criteria:  (a) 
what  information  is  required;  (b)  how  the  information  is  extracted  from 
measurements;  and  (c)  the  cost  and  time  involved  in  gathering  the 
information. 

2.  Test  programs  are  classified  as  either  destructive  or  nonde¬ 
structive,  depending  on  the  method  of  measurement.  In  a  destructive 
program,  samples  are  removed  at  specific  points  (usually  widely  spaced) 
and  are  tested  at  the  site  or  in  the  laboratory.  This  type  of  program 
can  be  very  quick  and  cheap  or  very  time-consuming  and  expensive,  de¬ 
pending  on  the  information  required  and  the  material  being  sampled. 

3.  Nondestructive  testing  allows  the  user  to  look  at  a  large 
portion  of  an  area,  usually  more  quickly  than  in  destructive  testing, 
and  does  not  require  the  removal  of  any  material.  However,  the  equip¬ 
ment  used  for  nondestructive  testing  is  quite  expensive,  and  the  data 
collected  are  difficult  to  interpret.  For  example,  seismic  surveys  are 
used  to  obtain  acceptably  accurate  measurements  of  layering  in  many 
situations,  but  detection  of  an  interface  between  two  layers  is  diffi¬ 
cult  or  impossible  if  a  high-propagation-velocity  material  overlies  a 
low -propagation -velocity  material.  Although  this  situation  does  not 
ordinarily  occur  in  natural  soils,  it  is  almost  invariably  present  in 
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artificial  media,  such  as  pavements. 

4.  The  above-mentioned  considerations  led  to  the  idea  that  mea¬ 
surements  made  with  microwave  energy  might  provide  an  acceptable  solu¬ 
tion  to  the  problem  of  nondestructive  testing  of  layered  media  such  as 
soils  and  rocks.  Microwave  energy  can  penetrate  most  terrain  materials 
to  some  degree  and  is  sensitive  to  changes  in  the  electrical  properties 
of  these  materials.  In  addition  to  being  nondestructive,  microwave  sys¬ 
tems  are  noncontact  devices,  a  fact  which  allows  a  wide  variety  of  mount¬ 
ing  options  ranging  from  air  platforms  to  surface  vehicles. 

5.  Deep  sensing  through  ice  with  radar  systems  operating  at  high 
frequencies  has  been  perfected.  Results  of  such  sensing  operations  show 
layers  within  the  ice  (i.e.  changes  in  ice  density  due  to  age)  as  well 
as  the  total  ice  depth.  Shallow  sensing  with  these  devices  Is  not  fea¬ 
sible  because  of  low  resolution;  also,  the  transmitted  signal  interferes 
with  the  near-surface  return  (i.e.  the  duration  of  the  transmitted  pul re 
is  too  long). 

6.  Sensing  through  soils  or  other  terrain  materials  demands  the 
use  of  low  frequencies  to  achieve  penetration.  In  general,  the  lower 
the  frequency,  the  less  attenuation  the  wave  experiences  when  traveling 
in  a  given  material.  The  duration  of  the  transmitted  pulse  must  be  kept 
short  to  provide  an  acceptable  degree  of  resolution  or  time  discrimina¬ 
tion  when  thicknesses  of  thin  layers  are  measured.  Systems  with  ex¬ 
tremely  short  pulses  (monopulse  systems)  operating  at  lov.T  frequencies 
have  not  been  used  very  successfully  because  of  data  interpretation 
problems  associated  with  antenna  propagation  (ringing  and  dispersion). 

7.  Swept -frequency  radar  can  provide  an  acceptable  depth  of  pen¬ 
etration  into  terrain  materials  when  the  transmitter  is  operating  at 
low  radar  frequencies.  Pulse  width  and  distortion  proolems  are  elimi¬ 
nated  by  operating  in  a  continuous  wave  mode.  Surface  and  subsurface 
information  can  then  be  extracted  through  interpretation  of  interference 
patterns  found  in  the  reflectance  data. 

Previous  WES  Investigations 

8.  Several  investigations  with  microwave  systems  have  already  been 
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conducted  at  the  U.  S.  Array  Engineer  Waterways  Experiment  Station  (WES) 

and  the  results  have  been  published.  These  investigations  included  a 

laboratory  study  to  examine  the  possible  use  of  pulsed  radar  systems 

to  detect  and  measure  the  depths  of  surface  water  and  groundwater  and 

1  2 

changes  in  soil  moisture  contents,  *  and  a  laboratory  study  with  a 

microwave  interferometer  that  measured  the  basic  electrical  properties 

of  various  test  samples  and  related  these  properties  to  the  physical 

o 

properties  of  the  test  sample.  The  results  of  these  studies  showed 
that:  (a)  the  electrical  properties  of  soils  and  other  materials  can 
be  measured  remotely  (the  range  used  in  the  radar  study  was  approximately 
lh  m*),  (b)  there  is  a  direct  correlation  between  a  soil’s  dielectric 
constant  and  water  content,  and  (c)  long -wavelength  microwave  energy 
will  penetrate  to  and  reflect  from  subsurface  interfaces. 

Purpose  and  Scope  of  Test  Program 

9.  The  purpose  of  the  overall  research  program  was  to  develop 
methods  for  extracting  quantitative  information  from  terrain  surfaces 
through  use  of  specially  designed  microwave  systems.  The  specific  ob¬ 
jective  of  this  study  was  to  determine  the  ability  of  a  swept -frequency 
radar  system  operating  under  field  conditions  to:  (a)  detect  the  pres¬ 
ence  of  substrata  of  various  thickness,  densities,  and  composition;  and 
(b)  develop  analytical  procedures  for  extracting,  from  the  reflectance 
curves,  information  related  to  the  electrical  properties  of  the  mate¬ 
rials  in  these  substrata. 

Definitions  of  Pertinent  Terms 


10.  Certain  terms  used  in  this  report  are  defined  as  follows: 
Attenuator.  A  device  that  reduces  the  amplitude  of  an  electrical 
signal  without  introducing  appreciable  phase  or  frequency  distortion. 


*  A  table  of  factors  for  converting  British  units  of  measurement  to 
metric  units  and  metric  units  of  measurement  .,0  British  units  is  pre¬ 
sented  on  page  ix. 
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Bandwidth.  The  difference  between  the  limiting  frequencies  of  a 
continuous  frequency  band.  The  bandwidth  of  a  device  is  the  difference 
between  limiting  frequencies  where  performance  falls  within  specified 
limits . 

Conductivity.  The  electrical  conductance  of  a  material  having 
unit  length  and  unit  cross  section. 

Dielectric  constant.  That  property  of  a  material  that  determines 
the  electrostatic  energy  stored  per  unit  volume  for  unit  potential  gra¬ 
dient;  synonymous  with  permittivity. 

Directional  coupler.  A  device  that  divides  the  input  signal  into 
two  output  signals.  Both  output  signals,  one  being  much  larger  than  the 
other,  are  porportional  to  the  input  signal. 

Incidence  angle.  The  angle  between  the  line  normal  to  the  surface 
at  the  point  of  incidence  of  the  electromagnetic  wave  and  the  line  of 
propagation  of  the  wave  to  the  surface. 

Isolator.  A  device  used  to  minimize  the  coupling  or  interaction 
of  components. 

L-band;  A  range  of  frequencies  from  0.39  to  1.5  GHz. 

Loss  tangent.  The  relation  between  the  dielectric  constant  and 
conductivity  at  a  given  frequency;  synonymous  with  dissipation  factor. 

Microwaves.  Radio  waves  that  have  wavelengths  so  short  that  they 
exhibit  some  of  the  properties  of  light. 

Octave.  The  interval  between  two  frequencies  having  a  ratio  of 
two  to  one. 

Optical  thickness.  The  distance  separating  two  reflecting  sur¬ 
faces  in  free  space  that  would  give  the  same  displacement  measurement 
as  that  obtained  from  the  two  surfaces  of  a  dielectric  layer. 

Polarization.  A  term  used  to  specify  the  direction  of  the  elec¬ 
tric  vector  in  a  polarized  radio  wave  with  respect  to  the  direction  of 
propagation. 

Power  reflection  coefficient.  The  ratio  of  the  reflected  power 
density  at  the  surface  of  a  material  to  the  incident  power  density. 
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PART  II:  TEST  PROGRAM 


Swept -Frequency  Radar  System 


11.  The  swept -frequency  radar  system  used  in  this  investigation 
(fig.  l)  was  designed  to  measure  reflectance-  from  a  test. site  over  the 
frequency  range  of  0.25  to  8.0  GHz.  The  system,  which  can  be  moved  from 
site  to  site  on  three 
vehicles,  is  composed 
of  two  major  sub- 
assemblies,  i.e.  the 
transmitter  and  the 
receiver  system. 

12.  The  trans¬ 
mitter  section  generates 
a  signal  directed  at  the 
test  site.  The  portion 
of  the  signal  reflected 
by  the  test  site  and  col¬ 
lected  at  the  receiver 
antenna  is  processed 
through  the  receiver's 
electronic  circuits  to  a 
recorder. 

13.  The  principal 
components  of  the  swept - 
frequency  radar  system 
are  shown  in  fig.  2.  De¬ 
tails  are  presented  in 
the  following  paragraphs 
on  how  the  various  func¬ 
tions  are  performed. 

The  heart  of  the  trans-  .  .  ,  ,  _ 

Fig.  1.  Truck-mounted  swept -frequency 

mitter  section  is  the  radar  system 
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Fig.  2.  Swept -frequency  radar  block  diagram 

0.25-  to  8.0-GHz  sweeping  oscillator.  On  command,  it  generates  fre¬ 
quencies  over  octave  bands  (sweeping  one  octave  at  a  time)  at  selected 
power  levels.  The  rated  power  output  for  each  octave  is  shown  in  the 
following  tabulation: 

Hated  Power  Output 

Frequency  Range,  GHz  milliwatts 
0.25  -  0.5  20 

0.5  -1.0  50 

1.0  -  2.0  60 

2.0  -  4.0  50 

4.0  -  8.0  25 


The  rate  of  frequency  change  can  be  controlled  by  means  of  an  external 
signal.  The  output  power  can  be  leveled  by  sampling  a  portion  of  the 
transmitted  power  at  a  control  point  and  routing  a  d-c  voltage 
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proportional  to  the  sampled  power  back  to  the  sampling  input  on  the 
sweeping  oscillator. 

14.  The  0.5-  to  S.O-GHz  signals  generated  by  the  sweeping  oscil¬ 
lator  are  routed  to  the  appropriate  traveling  wave  tube  amplifier  through 
an  isolator  (10-db  attenuator).  The  0.25-  to  0.5-GHz  signals  need  no 
further  amplification,  since  attenuation  through  the  various  components 
is  low  at  these  frequencies.  Traveling  wave  tube  amplifiers  are  used 
because  of  their  relatively  wide  bandwidth  and  uniform  gain  character¬ 
istics;  the  minimum  rated  gain  and  rated  power  output  are  given  below: 


Frequency 

GHz 

Minimum 

Rated 

Gain 

db 

Rated 

Power 

Output 

w 

0.5  -  1.0 

30 

1.0 

1.0  -  2.0 

20 

10.0 

2.0  -  4.0 

30 

10.0 

4.0  -  8.0 

30 

10.0 

Isolators  are  used  between  the  sweeping  oscillator  and  the  traveling 
wave  tube  amplifiers  to  adjust  the  power  levels  to  the  best  range  of 
operation  and  to  reduce  interaction  between  the  two  components. 

15.  The  signals  are  routed  from  the  traveling  wave  tube  amplifier 
(or  directly  from  the  sweeping  oscillator  in  the  case  of  the  0.25-  to 
0.5-GHz  octave)  through  the  single  transmitter  coaxial  cable  to  a  co¬ 
axial  switch.  The  coaxial  switch,  controlled  by  the  octave  selector, 
directs  the  signals  to  the  proper  power  sampler  (i.e.  one  of  five  di¬ 
rectional  couplers). 

1 6.  The  power  sampler  divides  the  input  power  into  two  components, 
one  much  larger  than  the  other.  The  smaller  component,  proportional  to 
the  signal  applied  to  the  transmitting  antennas,  is  rectified  with  a 
crystal  detector  and  routed  to  the  transmitter  level  amplifier.  Two 
outputs  are  obtained  from  this  amplifier;  one  is  available  for  record¬ 
ing  in  the  receiver  section,  and  the  other  is  fed  back  to  the  leveling 
input  on  the  sweeping  oscillator  to  keep  the  transmitted  power  at  a 
constant  level. 
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17.  The  larger  output  from  the  power  sampler  is  the  signal  ap¬ 
plied  to  the  transmitter  antenna.  A  coaxial  switch  controlled  by  the 
antenna  selector  connects  the  proper  power  sampler  to  the  0.25-  to 
4.0-GHz  transmitting  antenna.  The  4.0-  to  8.0-GHz  transmitting  antenna 
remains  connected  to  the  4.0-  to  8.0-GHz  power  sampler  at  all  times. 

18.  Two  parabolic  transmitting  antennas  are  used  to  propagate  the 
frequencies  generated  in  the  transmitter  section  through  spat  i  toward 
the  boresighted  center  of  the  top  surface  of  the  test  site.  The  two 
antennas  have  the  following  characteristics: 


Antenna  diameter,  m 

Antenna  polarization 

Antenna  feed 

Half -power  beam  width, 
deg,  at: 

0.3  GHz 

2.5  GHz 

5.0  GHz 

7.0  GHz 

Antenna  gain  (over  an  iso¬ 
tropic  radiator),  db,  at: 

0.3  GHz 

2.5  GHz 

5.0  GHz 

7.0  GHz 

19.  Two  parabolic  receiving  antennas  are  used  to  collect  a  portion 
of  the  microwave  energy  reflected  from  the  test  site.  They  are  identi¬ 
cal  with  the  two  transmitting  antennas  and  are  also  boresighted  toward 
the  center  of  the  top  surface  of  the  test  site. 

20.  Through  the  application  of  the  antenna  selector  switch,  the 
signal  from  the  proper  receiving  antenna  is  rectified  with  a  crystal  de¬ 
tector  and  amplified  in  the  receiver  amplifier.  The  amplifier  gain  is 
controlled  with  a  gain  switch  and  has  four  ranges:  0-30,  0-100,  0-300, 
and  0-1000  pw. 
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21.  The  receiver  signal  is  sent  to  the  impedance  matching  ampli¬ 
fier,  which  reduces  the  source  impedance  and  scales  the  amplitudes  for 
further  processing. 

22.  The  input-output  selector  is  used  to  select  either  the  re¬ 
ceiver  signal  from  the  impedance  matching  amplifier  or  the  transmitter 
signal  from  the  transmitter  level  amplifier.  The  selected  signal  is 
then  routed  either  to  an  x-y  recorder  for  a  graphic  display  or  to  the 
analog -to-digital  (A-D)  converter  for  further  processing.  An  oscillo¬ 
scope  is  also  used  to  monitor  the  input  signal  at  this  point. 

23.  The  A-D  converter  digitizes  the  input  signal  at  a  200- 
character-per-seeor.d  rate  on  command  from  the  interface  and  control  unit. 
The  digitizing  rate  is  well  within  the  capability  of  the  A-D  converter 
(bandwidth  in  excess  of  100,000  Hz)  and  is  fast  enough  to  stop  power 
fluctuations  caused  by  antenna  movement.  The  output,  a  12 -bit  binary 
number,  is  fed  to  the  interface  and  control  unit. 

24.  The  interface  and  control  unit  prepares  the  digital  signal 
for  recording  on  the  digital  tape  recorder.  The  recorder  has  two  modes 
of  operation: 

a.  A  fixed  data  mode,  in  which  data  from  12  thumb  wheel 
switches  can  be  recorded  on  the  tape. 

b.  An  external  data  mode,  in  which  the  data  from  the 
analog -to -digital -converter  can  be  recorded. 


The  fixed  data  mode  permits  identification  (site  number,  octave  fre¬ 
quency  band,  receiver  gain  setting,  transmitter  or  receiver  power  re¬ 
cording,  and  sample  or  reference  standard  recording)  of  the  external 
data.  In  the  external  data  mode,  256  data  words  are  recorded  for  each 
frequency  octave.  The  use  of  256  data  words  per  octave  was  adequate  to 
resolve  the  highest  frequencies  of  interest  in  the  data  record.  In  ad¬ 
dition,  each  time  a  digital  word  is  recorded,  a  signal  is  sent  from  the 
interface  and  control  unit  to  the  sweep  integrator. 

25.  The  sweep  integrator  conditions  the  pulse  from  the  interface 
and  control  unit  and  applies  the  resulting  signal  to  an  integrator  cir¬ 
cuit.  The  output  is  a  ramp  voltage  (an  analog  signal)  that  starts  at 


zero  volts  and  increases  an  equal  amount  each  time  the  signal  is  digi¬ 
tized  and  recorded.  This  ramp  voltage  is  fed  back  to  the  sweeping  os¬ 
cillator  for  frequency  control.  The  generated  frequency  is  lowest  when 
the  ramp  Voltage  is  zero  and  increases  to  twice  the  lowest  frequency 
(one  octave)  at  the  highest  ramp  voltage. 

Test  Sites 


2 6.  Three  sites  on  Interstate  Highway  20  near  Clinton,  Miss, 
were  selected  for  tests  with  the  swept-frequency  radar  system.  The 
highway  was  constructed  on  a  silt  foundation  and,  at  the  time  of 
the  tests,  had  an  8.5-in.  asphaltic  concrete  surface  (final  design 
thickness  of  the  asphaltic  concrete  was  11  in.).  The  pavement  was 
separated  from  the  silt  foundation  by  a  clay-gravel  base  course.  The 
profile  of  a  typical  highway  section  is  shown  in  fig.  3,  and  the  high¬ 
way  data  are  given  in  table  1.  A  summary  of  the  soils  data  collected  by 
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Fig.  3.  Highway  profile  for  westbound  traffic  (not  to  scale) 


the  Mississippi  State  Highway  Department  in  the  vicinity  of  the  test 
sites  is  given  below: 

a.  Clay-gravel  base.  This  granular  material  had  from  10  to 
26  percent  fines  (silt-  and  clay-sized  particles  passing 
a  Ro.  200  sieve)  ana  from  35  to  48  percent  gravel  (ma¬ 
terial  retained  on  a  No.  U  sieve).  It  classified  SP-SC 
.  to  Sid  according  to  the  Unified  Soil  Classification  System 

(USCS).  The  Atterberg  limits  for  material  passing  a  Ro. 

Uo  sieve  were:  liquid  limit,  18-23;  and  plasticity  index, 
4-9  (some  samples  were  nonplastic). 
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b.  Silt  foundation.  This  fine-grained  material  (99  percent 
”  fines)  at  site  1  (sta  724)  and  site  2  (sta  654)  clas¬ 
sified  CI>  according  to  the  USCS.  Its  average  Atterberg 
limits  were:  liquid  limit,  J3;  and  plasticity  index,  l4. 
At  site  3  (sta  6o4)  the  soil  had  less  clay  (96  percent 
fines)  and  classified  CL-ML.  Its  average  Atterberg  lim¬ 
its  were:  liquid  limit,  28;  and  plasticity  index,  6.  At 
site  1  the  finish  grade  was  the  natural  soil;  at  site  2 
approximately  1  m  of  fill  was  placed  over  a  1-m-high  by 
1.5-m-wide  concrete  box  culvert;  and  at  site  3  approxi¬ 
mately  0.6  m  of  fill  was  placed  to  attain  roadbed  eleva¬ 
tion.  The  fill  materials  were  the  same  as  the  silt 
foundation  except  for  being  reworked  and  compacted. 

27.  The  moisture  contents  and  dry  densities  of  the  foundation  ma¬ 
terials  varied  from  point  to  point  along  the  highway.  For  example,  at 
site  1  (sta  724),  during  highway  construction,  the  moisture  content  and 
dry  density  of  the  clay -gravel  base  were  4.4  percent  and  2.04  g/em  , 
respectively.  Moisture  content  in  nearby  sample  points  ranged  from 

O 

3.8  to  7  percent,  and  dry  density  from  1.94  to  2.17  g/em  .  The  varia¬ 
tions  were  similar  for  the  silt  foundation;  the  moisture  content  ranged 

O 

from  20  to  25  percent,  and  the  dry  densities  from  1.52  to  1.81  g/cnr . 
These  same  variation  ranges  also  were  noted  at  sites  2  and  3. 

28.  Because  of  the  above-mentioned  variations  in  density  and 
moisture  content,  a  single  relative  dielectric  constant  value  (neces¬ 
sary  for  later  analysis  of  data)  was  chosen  for  each  highway  material. 
These  dielectric  constants  were  representative  of  the  length  of  highway 
on  which  tests  were  made.  The  constants  were  obtained  from  curves  of 

relative  dielectric  constant  versus  water  content,  determined  from  pre- 

3 

vious  laboratory  investigations  at  the  VIES.  Two  such  curves,  shown  in 
fig.  4,  represent  data  obtained  from  414  samples  of  sands,  silts,  and 
clays  (actually  .12  different  soils). 

29.  For  the  clay -gravel  base,  a  moisture  content  of  5  percent  and 
a  dry  density  of  2  g/cnr  were  selected.  These  are  equivalent  to  a  water 
content  of  0.1  g/cnr  (the  product  of  moisture  content  and  dry  density), 
resulting  in  a  relative  dielectric  constant  of  6.5  (from  fig.  4).  For 
the  silt  foundation,  a  moisture  content  of  19  percent  and  a  dry  density 
of  1.6  g/cnr  were  selected.  These  are  equivalent  to  a  water  content  of 
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0.3  g/cm  ,  resulting  in  a  relative  dielectric  constant  of  approximately 
1 6.  The  electrical  properties  of  the  asphaltic  concrete  were  obtained 
from  laboratory  measurements  with  an  L-band  interferometer  on  specially 
prepared  84-cm-square  by  5-cm-deep  samples  of  coarse  asphaltic  concrete 
(coarse  aggregate  that  can  pass  through  a  l-l/2-in.  screen).  This  sam¬ 
pled  material  was  identical  with  the  asphaltic  concrete  at  the  highway 
sites;  for  the  laboratory  tests,  the  angle  of  incidence  was  k5  deg  (fixed 
by  interferometer  design)  and  the  temperature  was  2k  C,  or  room  tempera¬ 
ture.  A  relative  dielectric  constant  of  3 • 53  was  selected  for  this  mate¬ 
rial.  The  measured  electrical  properties  are  given  in  the  tabulation 

below :  n  ,  . .  _  , 

Relative  Conduc- 


Frequency 

Dielectric 

Loss 

tivity 

GHz 

Polarization 

Constant 

Tangent 

mho/m 

1.07k 

Parallel 

3.55 

0.02k 

0.005 

Perpendicular 

3.55 

0.038 

0.008 

1.12k 

Parallel 

3.k9 

0.013 

0.003 

Perpendicular 

3.58 

0.032 

0.007 

1.311 

Parallel 

3.k7 

0.021 

0.005 

°erDendicular 

3.52 

0. 0k9 

0.012 

1.102 

Parallel 

3.51 

0.009 

0.002 

Perpendicular 

3.52 

0.0k8 

0.008 

l.k9? 

Parallel 

3.5k 

0.013 

o.ook 

Perpendicular 

3.52 

0.025 

0.007 

Reference  Standard 

30.  The  reflected  signal  from  a  reference  standard  was  measured 
prior  to  each  highway  measurement.  Any  changes  in  system  operation 
showed  up  as  changes  in  both  measurements  and  were  canceled  when  the 
power  reflectance  was  calculated  (see  paragraph  35) • 

31.  The  reference  standard  was  constructed  of  metal  plates  sup¬ 
ported  on  a  metal  frame  in  a  fold-up  design  (fig.  5)  that  made  the  stan¬ 
dard  portable.  Because  the  high -conductivity  metal  made  a  microwave 
mirror,  nearly  all  the  incident  power  was  reflected  (or  reradiated) 
specularly. 

32.  When  set  up,  the  reference  standard  measured  3.35  m  square. 
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Fig.  5.  Reference  standard 


There  was  no  middle  seam  in  the  metal  plate  to  interfere  with  reflectance 
measurements  made  when  the  antenna  beam  widths  were  stall.  When  the  an¬ 
tenna  beam  widths  were  large,  there  were  four  seams  in  the  Illuminated 
area.  However,  since  the  seam  area  was  small  compared  with  the  total 
illuminated  area,  there  was  very  little  effect  on  the  total  reflected 
signal. 

Test  Procedures 


33.  It  was  important  that  each  swept -frequency  radar  test  be  per¬ 
formed  in  an  identical  fashion  so  that  variations  in  the  data  would  be 
caused  by  the  reflectance  variations  of  the  test  site  alone.  Procedures 
were  developed  to  '  :sure  that  the  relative  positions  of  the  antennas 
with  respect  to  the  reference  stanu^rd  (when  in  place)  would  remain  the 
same  even  though  the  equipment  was  moved  from  site  to  site.  These  pro¬ 
cedures  are  described  in  the  following  paragraphs. 


34.  The  midpoint  between  the  transmitting  and  receiving  antennas 
was  positioned  directly  above  the  test  site.  The  antennas  were  rotated 
with  electric  motors  about  the  midpoint  until  they  were  in  a  horizontal 
plane,  and  the  plane  formed  by  the  lines  of  bores ight  was  vertical.  This 
leveling  procedure  was  checked  with  level  lights  controlled  by  mercury 
switches  that  turned  on  when  the  antennas  were  in  prop*..'  alignment.  The 
range  (the  vertical  height  of  the  antenna  carriage  minu.--  bhe  distance 
the  antennas  were  mounted  below  the  carriage)  used  for  these  tests  was 
7.6  m.  The  reference  standard  was  moved  to  the  center  of  the  test  site 
and  leveled  by  using  a  bubble  level  at  each  corner.  The  position  of  the 
reference  standard  in  the  center  of  the  test  site  was  checked  with  a 
closed-circuit  television  system,  the  camera  of  which  was  mounted  at  the 
midpoint  between  the  transmitting  and  receiving  antennas  on  the  antenna 
carriage  in  boresight  with  the  antennas.  The  position  and  height  of  the 
antennas  were  also  checked  with  a  weighted  line  connected  at  the  same 
point  as  the  television  camera  midway  between  the  transmitting  and  re¬ 
ceiving  antennas. 

35.  The  transmitter  and  receiver  sections  of  the  system  were  en¬ 
ergized,  and  the  transmitted  power  and  the  power  reflected  by  the  ref¬ 
erence  standard  to  the  receiving  antenna  as  the  frequency  was  varied 
were  recorded.  Although  the  power  applied  to  the  antennas  was  main¬ 
tained  at  a  constant  value  over  each  octave  band  (see  paragraph  16), 
the  antennas  were  not  uniformly  efficient  at  all  frequencies,  and  the 
power  collected  by  the  receiving  antennas  from  the  reference  standard 
varied  in  amplitude  (fig.  6).  The  reference  standard  was  then  removed, 
exposing  the  surface  of  the  test  site.  The  transmitted  and  received 
power  were  again  recorded.  The  sample  power  reflectance  referenced  to 
the  standard  reflectance  was  calculated  for  all  frequencies  as: 

P  /P  P  P 

D  _  rRS/rTS  _  rRS  TR 

K  P  /P  p  p  ”  ' 

RR/  TR  RR  TS 

R  =  power  reflection  coefficient 
P_„  =  reflected  power  froir,  che  test  site,  w 


where 


Prpg  =  transmitted  power  with  the  test  site  exposed,  w 
=  reflected  power  from  the  referenced  standard,  w 

PTR  =  transmitted  power  with  the  reference  standard  in  place,  w 

36.  Three  replicate  tests  were  made  at  each  test  site.  The  data, 
in  terms  of  power  reflectance,  were  then  used  for  layer  analysis  and  for 
the  calculation  of  electrical  properties  (see  paragraphs  6O-69). 

Test  Results 

37.  The  curves  of  power  reflectance  versus  frequency  are  given  in 
plates  1-3  for  site  1,  plates  4-6  for  site  2,  and  plates  7-9  for  site  3* 
Each  of  these  curves  was  developed  by  applying  equation  1  to  the  data 
obtained  from  the  four  separate  recordings  of  transmitted  and  received 
power  for  both  the  reference  standard  and  the  test  site. 

38.  The  accuracy  of  the  measured  data  must  be  considered  to  be 
better  at  some  frequencies  than  at  others.  For  example,  the  beam  of 
the  transmitting  and  receiving  antennas  at  the  lowest  operating  fre¬ 
quencies  (around  0.25  GHz)  was  wide  enough  at  a  height  of  7.6  m  to  allow 
part  of  the  projected  beam  to  fall  outside  the  reference  standard.  Also, 
the  antennas  used  in  this  system  had  limited  bandwidths.  The  low- 
frequency  antennas  had  a  specified  bandwidth  of  0.3  to  3.0  GHz,  and  the 
high-frequency  antennas  had  a  bandwidth  of  approximately  4.5  to  8.0  GHz. 
Microwave  energy  outside  these  frequency  ranges  was  not  handled  effi¬ 
ciently,  as  shown  by  the  extremely  low  amplitudes  recorded  for  these 
points  in  fig.  C.  Arithmetic  operations  on  such  data  tend  to  have  a 
higher  percentage  of  error  associated  with  each  calculation  than  those 
on  data  inside  the  antenna  bandwidths.  For  these  reasons,  analysis  was 
concentrated  on  data  in  the  0.5-  to  2.0-GHz  frequency  region.  The 

4.5-  bo  S.O-GHz  frequency  region  provided  information  on  surface  and 
near-surface  properties  because  of  its  lesser  penetrating  power.  The 
0.25-  to  0.5-GHz  frequency  region  had  the  highest  penetrating  power, 
but  v:as  used  only  to  check  for  trends  because  of  the  errors  discussed 
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PART  III:  REDUCTION  AND  ANALYSIS  OF  DATA 

Target  Model 

39*  When  an  electromagnetic  wave  traveling  in  one  medium  impinges 
I  upon  a  second  medium  having  differing  electrical  properties,  the  wave  gen¬ 

erally  will  be  partially  transmitted  and  partially  reflected.  The  trans- 
;  mitted  wave  is  attenuated  as  it  travels  through  the  second  medium  and, 

upon  encountering  a  third  medium,  experiences  partial  reflection  and 
transmission  again.  Theoretically,  this  process  can  continue  indefinitely. 
The  reflected  signals  of  interest  are  those  that  travel  back  to  the  sur- 
i  face  with  enough  amplitude  to  significantly  alter  the  total  amplitude  of 

the  return  signal.  The  amplitudes  of  the  reflected  signals  are  governed 
by  the  properties  of  the  medium.  Certain  assumptions  were  made  regarding 
these  properties  in  order  to  simplify  the  methods  used  to  extract  infor¬ 
mation  about  the  medium  from  the  reflection  curves,  as  discussed  below. 

40.  The  test  sites  were  assumed  to  be  composed  of  layers  of  var- 
i  ious  thicknesses.  The  material  in  each  layer  was  uniform  and  did  not 

have  any  discontinuities  either  within  the  layer  or  at  its  top  and  fcot- 
S  tom  surfaces.  The  electrical  properties  of  each  material  were  assumed 

to  be  frequency  independent  (at  least  over  the  frequency  range  used  for 
calculations)  and  predictable  from  past  tests  in  the  laboratory. 

|  41.  The  transition  zones  between  the  layers  were  assumed  to  be 

very  thin  (negligible)  and,  thus,  the  layer  interfaces  abrupt.  A  thick 
|  transition  zone  would  tend  to  reduce  the  amplitude  of  the  boundary  re- 

I  flection  and  might  cause  a  layer  to  be  overlooked.  This  is  a  matter  of 

|  degree,  since  the  thickness  is  measured  relative  to  the  wavelength  of 

i  the  transmitted  signal  in  the  material  through  which  the  wave  is  travel¬ 

ing.  A  thin  transition  zone  at  long  wavelengths  might  be  a' very  thick 
transition  zone  at  short  wavelengths.  For  the  purpose  of  this  report,  a 
thin  transition  zone  is  one  whose  thickness  is  less  than  one-eighth  of 
’  a  wavelength  of  the  transmitted  signal  (Rayleigh  criterion). 

42.  Each  layer  surface  was  assumed  to  be  smooth.  Again,  smooth¬ 
ness  is  a  matter  of  degree  and  can  be  evaluated  by  comparing  the  surface 
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roughness  deviations  to  one -eighth  of  a  wavelength  of  the  transmitted 
signal  in  the  same  manner  as  that  for  the  transition  zones.  The  assump¬ 
tion  for  smoothness  allows  prediction  of  surface  electrical  properties, 
since,  in  general,  roughness  tends  to  reduce  the  signal  level  at  the 
radar  receiver  by  scattering  the  reflected  energy  and  thus  tends  to  de¬ 
crease  the  accuracy  of  the  computations. 

43.  In  addition,  the  top  and  bottom  surfaces  of  each  layer  were 
assumed  to  be  parallel  and  perpendicular  to  the  radar  beam.  If  the 
layers  were  not  parallel,  the  layer  thickness  would  not  be  uniform. 

Also,  a  sloped  surface  would  tend  to  direct  the  reflected  energy  away 
from  the  receiver  and  thus  reduce  the  amplitude  of  that  signal  component 

44.  These  assumptions  do  not  restrict  operation  to  only  ideal 
test  sites;  rather,  the  accuracy  of  computed  results  is  increased  when 
the  assumptions  are  met. 


Spectral  Response 


45.  The  reflectance  of  a  material  can  be  highly  dependent  on  the 
frequency  of  the  microwave  source.  Materials  are  roughly  divided  into 
two  classes:  (a)  conductors  and  (o)  dielectrics,  or  insulators.  The 
dividing  line  between  the  two  classes  is  not  sharp,  and  some  materials 
(soil,  for  example)  are  considered  to  be  conductors  in  one  part  of  the 
radio  frequency  range  and  dielectrics  in  another  part.  The  criterion 
for  classifications  is  the  value  of 


a 

where 

c  =  electrical  conductivity,  mhos  m 
<o=  circular  frequency  (  =  2nf),  radians/'sec 
f  =  frequency,  Hz 
e  =  dielectric  constant,  farads/m 

When  ct4o€  (which  is  also  defined  as  the  ratio  of  conduction  current 
density  to  displacement  current  density  in  the  material)  is  greater 
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than  1,  the  material  is  a  conductor,  and  when  a /me  is  less  them  1,  the 
material  is  a  dielectric.  At  high  frequencies  (i.e.  above  a  few  mega¬ 
hertz),  most  nonmetallic  terrain  materials  act  as  dielectrics  insofar 
as  their  reflecting  properties  are  concerned.  The  waves  transmitted 
through  the  surface  into  the  material  are,  however,  rapidly  attenuated. 

1+6.  The  measurements  discussed  in  this  report  were  made  at  high 
frequencies  (0.25  to  8.0  GHz)  and,  therefore,  the  materials  are  con¬ 
sidered  to  be  dielectrics,  or  poor  conductors.  For  these  materials 
(especially  low-loss  materials),  energy  penetration  is  possible,  and  re¬ 
flection  from  subsurface  layers  can  combine  with  surface  reflections  tc 
give  a  total  amplitude  highly  dependent  on  the  layer  geometry  and  trans¬ 
mitter  frequency.  Thus,  a  measurement  at  a  single  frequency  cannot  be 
used  to  predict  the  surface  reflectance  unless  the  test  site  is  known 
to  be  without  layers. 

1+7.  In  general,  the  reflectance  from  an  unlayered  test  site  will 
change  with  increasing  frequency  because  of  the  decreasing  conductivity 
effect.  Thus,  in  some  cases,  electrical  properties  can  be  predicted  by 
making  measurements  at  two  frequencies,  one  much  lower  than  the  other, 
and  noting  the  change  in  measurement  amplitudes.  This  is  not  possible 
for  most  field  measurements  since  test  sites  usually  have  some  degree 
of  layering. 

1+8.  For  layered  test  sites,  subsurface  reflection  can  combine 
with  the  surface  reflection  to  produce  a  per-.odic  amplitude  change  as 
the  frequency  is  varied.  This  amplitude  change  is  caused  by  destruc¬ 
tive  and  constructive  wave  interference  and  decays  with  increasing 
depth.  In  the  case  of  field  measurements,  some  scattering  occurs  from 
rough  surfaces  and  in  the  transition  zones  between  subsurface  layers . 
This  scattering  reduces  the  influence  of  multiple -reflected  waves.  The 
net  result,  an  increase  in  the  noise  level  of  the  reflected  signal, 
gives  a  somewhat  random  appearance  to  the  reflectance  curves . 

Spectral  Analysis 

1+9.  The  following  analysis  outlines  the  derivation  of  the 
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a  =  attenuation  constant,  m”1 
P  =  phase  constant, 

The  propagation  factor  can  he  computed  from 

7  =  +  jme) 


'  (3) 


where 

H  =  magnetic  permeability  (p.  =  HrU0)>  H/m. 

e  =  dielectric  constant  (e  =  er€0),  farads/m 

(j  =  relative  magnetic  permeability,  dimensionless 
r  rt 

4  =  free  space  magnetic  permeability,  4tt  X  10  '  H/m 
=  relative  dielectric  constant,  dimensionless 
€q  =  free  space  dielectric  constant,  1/(36tt  X  10^),  farads/m 
The  attenuation  constant  a  is  equal  to  the  real  part  of  the  propaga¬ 
tion  factor: 


a  =  co 
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The  phase  constant  3  is  equal  to  the  imaginary  part  of  the  propaga¬ 
tion  factor: 


(5) 


For  wave  propagation  in  low-loss  and  nonmagnetic  media 

(a  c  V  (a  e 

and  rr- 


(l  +  ^a) 
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c  V  r 
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where 


c  =  free-space  wave  velocity,  300  x  10  m/sec. 


23 


50.  If  the  frequency  is  sufficiently  high,  P  »  a  .  The  power 
reflections  can  be  simplified  from  the  more  general  expression  found  in 
the  literature  (see,  for  example,  reference  4)  to  the  following : 


Ei“ 


VS 


cos  0 
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er2  -  £rl  sin2  9 , 


1/2 


Vs i cos  6  +  (6 


r2  “  erl  sin2  8 


1 72 
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where  subscripts  1  and  2  on  refer  to  media  1  and  2  in  a  layered 
system  and 

R j^=  power  reflection  coefficient  for  perpendicular  polarization 
R^  =  power  reflection  coefficient  for  parallel  polarization 
0  =  angle  of  incidence,  deg 


51.  The  reflection  coefficients  versus  incidence  angle  calculated 
from  equations  8  and  9  for  fresh  water  (*r2  =  8l)  and  dry  soil  (ey2  =  6) 
with  medium  1  being  air  =  1)  are  shown  in  fig.  8.  As  has  been 
stated  (paragraph  45),  at  frequencies  greater  than  a  few  megahertz, 
these  media  act  as  dielectrics  (low-loss  media)  insofar  as  their  re¬ 


flecting  properties  are  concerned.  The  transmitted  wave  is,  however, 
rapidly  attenuated  by  the  effects  of  conductivity  (i.e.  in  the  factor  a). 
When  0=0 
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and  when  medium  1  is  air  .  =  l) 
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Fig.  8.  Reflection  coefficients  of  plane  radio  waves  from  soil 
and  freshwater  surfaces  at  frequencies  greater  than  10^  Hz 

52.  Note  that  the  reflectance”  in  fig.  8  does  not  change  signifi¬ 
cantly  with  incidence  angle  until  fl  reaches  10  or  15  deg;  therefore, 
the  slight  increase  in  incidence  angle  caused  by  the  separation  of  the 
transmitting  and  receiving  antennas  caused  only  minor  changes  in  the 
measured  reflectance  values. 

53.  For  a  layered  medium,  subsurface  reflections  can  combine  with 
surface  reflection  and  give  widely  varying  results,  depending  on  the 
thicknesses  of  the  layers,  the  layers’  electrical  properties,  and  the 
wavelength  of  the  transmitted  signal.  Consider  the  case  of  plane  waves 


at  normal  incidence  to  a  target  consisting  of  three  media  with  three 
propagation  factors  (fig.  9)* 


e, 


MEDIUM  1 
x  =  0 


MEDIUM  2 
x  =  d 


MEDIUM  3 


Fig.  9.  Reflection  and  transmission  of  plane  waves 
by  a  plane  sheet  at  normal  incidence 


54.  The  electrical  wave  reflectance  r  in  medium  1  can  be 
5  0 

written  as 
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(12) 


where  Er  and  axe  reflected  and  incident  electric  field  intensity 
respectively;  r^  and  Eire  electric  wave  reflections  at  the 

boundaries  between  media  1  and  2,  and  2  and  3,  respectively,  as  com¬ 
puted  from 


and  d  is  layer  thickness.  Impedances  Z  are  computed  for  each  layer 
as 


55.  When  equation  12  is  multiplied  by  its  complex  conjugate,  it 
becomes  an  equation  for  power  reflectance. 
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where 


2 

r,? |  =  power  reflection  coefficient  at  boundary  between 

media  1  and  2 


r2_ |  =  power  reflection  coefficient  at  boundary  between 

media  2  and  3 


and  9^  and  e23  are  defined  by 
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56.  As  shown  in  a  plot  of  equation  16  versus  layer  thickness  mea¬ 
sured  in  wavelengths  (fig.  10),  the  power  reflectance  can  have  a  number 


Fig.  10.  Power  reflectance  o 
at  normal  incidence  as  a  £ 
function  of  layer  thickness  k 
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of  values  at  one  single  frequency,  depending  on  the  layer  thickness  d  . 

57.  The  power  reflectance  also  varies  with  frequency  when  the 
layer  thickness  remains  constant  (fig.  11a).  The  positions  of  the 
maxima  and  minima  of  equation  1 6  for  low-loss  materials  can  he  found  by 
differentiating  the  reflectance  expression  with  respect  to  frequency 
and  by  setting  the  results  equal  to  zero  (the  value  for  slope  at  a  max¬ 
imum  or  minimum).  The  following  results  are  thus  obtained: 
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Thus,  a  maximum  or  minimum  occurs  for  the  power  reflectance  (equation  1 6) 
when 


sin  (20gd  -  023)  cos  0^  (l  -  R^) 
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This  expression  is  shown  graphically  in  fig.  11b  where  solution  points 
occur  at  the  intersections  of  the  two  curves.  Although  the  distance 
from  maximum  to  minimum  (see  fig.  ll)  is  different  from  the  distance 
from  minimum  to  maximum,  the  period  T^,  (the  distance  between  either 
adjacent  maxima  or  minima)  is  single-valued  and  is  displayed  when  20d 
increases  by  2tt  .  From  the  .results  of  equation  7 
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and 
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58.  From  equation  19,  the  layer  thickness  d  can  he  obtained 

from  the  cyclic  period  Tf  and  the  wave  velocity  of  that  layer. 

Layer  thickness  can  be  estimated  even  when  the  minima  and  maxima  are 
shifted  slightly  from  their  proper  frequency  positions  because  only  the 
period  of  oscillation  is  required.  Position  shifts  can  be  caused  by 
rough  surfaces,  transition  zones  between  materials,  nonhomogeneous  mate¬ 
rials,  effects  of  conductivity,  and  frequency  dependence  of  the  mate¬ 
rial's  electrical  properties.  The  accuracy  of  Tf  is  affected  by  the 
same  things  causing  position  shifts.  However,  wave  interference  phe¬ 
nomena  degrade  under  these  conditions,  and  before  a  significant  reduc¬ 
tion  in  Tf  is  observed,  the  cyclic  pattern  is  no  longer  identifiable. 

59.  The  results  of  equation  19  can  be  easily  extended  to  multiple- 
layer  sites  where  reflection  from  each  interface  causes  characteristic 
periodic  oscillations  to  appear  in  the  total  return.  These  periodic  os¬ 
cillations  can  be  extracted  by  means  of  a  correlation  process.  Repeated 
application  of  equation  19  can  be  used  to  determine  these  layer  thick¬ 
nesses  when  optical  depth  is  used  to  replace  the  physical  depth.  The 
optical  depth  can  be  defined  as  the  product  of  the  physical  depth  and 
the  square  root  of  the  relative  dielectric  constant.  Calculations  of 
these  types  were  employed  in  the  discussion  presented  in  paragraphs  63-70 . 
If  the  surface  and  subsurface  components  are  to  be  separated  from  the 
total  return,  equation  16  must  be  carried  one  step  further.  An  approxi¬ 
mate  expression  for  the  maxima  of  equation  16  can  be  found  by  setting 

cos  20d  »  1  (see  fig.  11 )  as  shown  below. 
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The  approximate  expression  for  the  minima  of  equation  lb  is  found  in  a 
similar  manner  by  setting  cos  2£d  =  -1  : 
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Then, 
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Note  that  if  the  subsurface  component  e  is  missing  or  neg¬ 
ligible,  R  =  R  .  -  R  ,  and  the  surface  component  is  R  =  R,0 
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Discussion  of  Results 


Surface  reflectance 

60.  The  level  of  the  power  reflectance  curve  is  established  by  the 
surface  reflectance  from  a  homogeneous  target  material,  or  by  the  largest 
signal  component  from  a  layered  site.  If  the  surface  reflectance  has  a 
larger  amplitude  than  any  of  the  subsurface  components,  the  total  sub¬ 
surface  reflectance  will  cycle  around  the  surface  reflectance.  If  a 
subsurface  component  has  the  largest  amplitude,  the  surface  reflectance 
will  cycle  around  the  subsurface  component  reflectance. 

61.  For  the  sites  measured  for  this  report,  the  surface  reflec¬ 
tance  was  much  larger  than  any  of  the  subsurface  components  and  was 
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assumed  to  be  the  averpge  power  reflectance  over  the  range  of  0.5  to 
2.0  GHz  (an  approximation  based  on  equation  24).  This  value  was  then 
substituted  into  equation  10  to  calculate  the  relative  dielectric  con¬ 
stant  of  the  surface.  The  results  of  these  calculations  are  shown  in 
table  2.  The  relative  dielectric  constants  varied  from  2.89  to  3*38 
for  the  nine  tests  and  averaged  3.12.  This  compares  very  well  with  the 
average  value  of  3.53  for  the  relative  dielectric  constant  from  L-band 
interferometer  measurements  (see  paragraph  29).  Note  that  the  highway 
surface,  sloped  away  from  the  traffic  center  line  (see  fig.  3).  The 
sloping  surface  of  the  test  site  would  have  the  same  effect  as  would 
operating  at  incidence  angles  other  than  0  deg  (see  fig.  8).  This  ef¬ 
fect  was  ignored  in  the  analysis  of  data  but  could  be  reduced  by  oper¬ 
ating  with  an  antenna  with  a  smaller  beam  width  or  by  operating  at  a 
lower  antenna  height  so  that  the  illuminated  area  on  the  pavement  sur¬ 
face  was  at  the  level  point  between  the  two  traffic  lanes. 

62.  The  antenna  height,  however ,  cannot  be  reduced  indefinitely. 
For  a  given  antenna  operating  at  a  given  frequency,  a  point  will  be 
reached,  as  the  height  is  reduced,  where  the  power  density  on  the  tar¬ 
get  area  is  not  uniform  or  easily  predictable  (i.e.  the  near  field  of 
the  antenna).  Although  some  of  the  data  presented  in  this  report  were 
actually  measured  in  the  near  field  of  the  antenna,  accuracy  was  not 
sacrificed,  since  the  refle  stance  of  the  test  site  was  referenced  to 

the  standard  reflectance,  with  both  at  the  same  range  (see  paragraph  35). 
Also,  since  the  layer  thicknesses  were  small  compared  with  the  antenna 
range  (or  height),  plane  w ave  analysis  could  still  be  used. 

Subsurface  reflectances 

63 .  Cyclic  patterns  are  present  in  all  the  power  reflectance 
curves  (plates  1-9).  Curves  obtained  from  measurements  at  the  same 
site  have  \;he  same  general  appearance  over  the  frequency  range  of 
0.25  to  2.0  GHz  and  differ  only  slightly  from  site  to  site.  Tne  curves 
have  high-frequency  signal  components  (noise)  riding  on  top  of  the  low- 
frequency  cyclic  patterns.  The  noise  does  not  appear  to  contain  in¬ 
formation  about  the  test  site  and  is  nrobably  the  result  of  digitizing 
procedures  and  minor  changes  in  transmitted  power. 


64.  In  order  to  extract  layer  thickness  from  the  reflectance  data, 
the  period  of  the  oscillations  had  to  be  determined.  This  was  done  by 
separating  the  original  power  reflectance  curve  into  its  various  signal 
components.  This  procedure  is  illustrated  in  fig.  12,  where  two  peri¬ 
odic  waves  of  decreasing  amplitudes  and  a  noise  component  are  extracted 
from  the  original  power  reflectance  curve.  The  periods  of  the  waves 
used  in  this  reconstruction  process  can  be  used  in  equation  19  to  com¬ 
pute  layer  thicknesses. 

65.  Another  method  used  to  extract  the  period  of  oscillation  in 
addition  to  the  graphical  technique  described  above  was  the  cross¬ 
correlation  technique.  The  data  were  cross -correlated  with  a  sine  wave, 
and  the  correlation  coefficient  showed  an  increased  amplitude  when  the 
sine  wave  period  was  the  same  as  that  of  the  periodic  data.  A  fast 
Fourier  transform  (?FT)  computer  program  with  the  following  modifica¬ 
tions  was  used  to  compute  the  correlations  (see  Appendix  A): 

a.  The  average  of  the  power  reflectance  data  (384  points 
equally  spaced  over  0.5  to  2.0  GHz)  was  computed  and 
subtracted  from  each  data  point. 

b.  Data  points  with  zero  amplitude  were  added  to  raise  the 
total  number  of  input  data  points  to  2048  (FFT  program 
requires  that  the  number  of  data  points  be  a  binary  num¬ 
ber  such  as  512,  1024,  ?.0k8,  etc).  The  number  2048  was 
selected  so  that  the  cross -correlation  plots  would  have 
a  smooth  continuous  appearance  when  plotted  by  the  com¬ 
puter.  This  gives  an  equivalent  frequency  range  to  the 
data  record  of  8.0  GHz. 

c_.  The  amplitudes  of  the  ITT  coefficients  were  multiplied  by 
5.333  (the  ratio  of  2048  to  384)  to  correct  for  the  ef¬ 
fect  of  the  added  zero -amplitude  points. 

d.  The  first  five  harmonics  were  not  used  (i.e.  the  first 
five  harmonics  had  periods  longer  than  the  data  record 
and,  therefore,  would  not  have  computed  properly). 

•'■6m  The  results  of  this  correlation  process  are  shown  in 
plates  10-12,  for  which  the  optical  depth  was  computed  as  follows: 

Optical  depth  in  meters  =  =  0.0l875n  (25) 

where  n  is  the  harmonic  number  from  the  FFT  computation.  Since  there 
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were  two  layers  (and,  therefore,  two  transition  zones)  in  the  highway- 
foundation,  each  plot  of  correlation  coefficient  versus  optical  depth 
showed  two  peaks.  The  first  peak,  at  approximately  0.1*  m,  was  easily 
identified  in  each  plate  and  indicated  the  optical  depth  to  the  bottom 
of  the  asphaltic  concrete  layer  (table  2).  The  second  main  peak  indi¬ 
cated  the  optical  depth  to  the  second  interface  (i.e.  between  the  clay- 
gravel  base  material  and  the  silt  foundation).  This  peak  was  easily 
identified  as  the  second  largest  peak  in  plate  10  for  site  1.  The  sec¬ 
ond  peak  for  site  2  (plate  11 )  was  also  taken  as  the  second  largest  peak, 
although  there  appeared  to  be  a  double  peak  associated  with  this  inter¬ 
face.  The  second  peak  for  site  3  (plate  12)  was  somewhat  more  difficult 
to  pick  out  with  only  a  single  measurement  being  used;  however,  the  lo¬ 
cation  for  the  second  peak  was  found  based  on  shape  similarity  when  all 
three  curves  were  compared. 

67.  The  remaining  peaks  present  in  plate  10  were  probably  caused 
by  the  analysis  process  itself;  short  data  records  can  produce  small 
peaks  in  the  correlation  curve  at  certain  frequencies,  depending  on  the 
character  of  the  data.  Long  data  records  without  sine  wave  distortion 
reduce  the  amplitudes  of  the  small,  ambiguous  correlation  peaks. 

68.  The  averages  of  the  correlation  curves  from  each  site  over 
the  optical  depth  range  of  0-2  m  are  shown  in  fig.  13.  The  amplitudes 
of  the  small,  ambiguous  correlation  peaks  were  reduced  from  those  of  the 
individual  curves,  since  they  were  displaced  fro;.*  curve  to  curve.  The 
peaks  that  were  due  to  interface  reflectances  were  maintained,  since 
they  were  found  at  nearly  the  same  optical  depth  in  each  measurement. 

The  layer  thicknesses  were  calculated  by  use  of  the  following  procedures: 

a.  Determine  the  optical  thickness  of  the  asphaltic  concrete 
layer  (i.e.  the  optical  depth  location  of  the  maximum  of 
•the  first  peak). 

b.  Calculate  the  asphaltic  concrete  thickness  by  dividing 
the  optical  thickness  by  the  square  root  of  the  asphal¬ 
tic  concrete  relative  dielectric  constant  (e  =3.53). 

r 

c.  Determine  the  optical  thickness  of  the  clay -gravel  base 
layer  (i.e.  the  difference  in  optical  depth  location 
between  the  maxima  of  the  first  and  second  peaks). 

d.  Calculate  the  clay-gravel  thickness  by  dividing  the 
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optical  thickness  by  phe  square  root  of  the  clay-gravel 

relative  dielectric  constant  (e  =6.5). 

r 

69.  The  results  of  these  proce<lures  are  tabulated  below  for  the 
average  data  in  fig.  13.  The  values  for  estimated  thickness  from  ground 
measurements  are  taken  from  bable  1. 

_ Swept  -Frequency^ ife,d.ar  Measurements _ 

Site  Optical  Predicted  Thickness  Estimated  from 
No.  Thickness,  m  Thickness,  in  Ground  Measurements,  m 


Asphaltic_  Concrete 


1 

0.455 

0.24 2 

0.226 

2 

0.415 

0.221 

0.210 

3 

0.380 

0.202 

0.214 

Clay-Grave^  pase 

1 

0.465 

0.182 

0.133  -  0.186 

2 

0.670 

0.263 

0.209  -  0.262 

3 

0.790 

0.310 

0.309  “  0.362 

Note  that  for  all  three  sites,  the  difference  between  the  radar  measure¬ 
ments  and  ground  measurements  for  the  Asphaltic  concrete  layer  was 
1.6  cm  or  less  (i.e.  e::*or  of  7  percePt  or  less).  For  the  clay-gravel 
base ,  the  radar  measurement  fell  withlh  t"he  range  of  ground  measurements 
for  sites  1  and  3  and  was  outside  the  fraPge  of  values  for  site  2  by  only 
1  mm. 

70.  If  the  same  techniques  were  uSed  to  predict  thickness  of 
layers  from  each  of  the  single  reflectance  tests  (as  listed  in  table  2), 
some  decrease  in  the  accuracy  of  the  computed  thicknesses  would  occur, 
as  shown  in  the  variation  from  test  t°  test  at  each  site.  For  example, 
predicted  thicknesses  of  the  asphaltic  concrete  varied  from  O.239  to 
0.2l \2  m,  0.221  to  0.226  m,  and  0.199  to  0.210  m  for  sites  1,  2,  and  3, 
respectively.  Predicted  thicknesses  ot  bhe  clay-gravel  base  varied  from 
0.182  to  0.188  m,  0.220  to  0.276  m,  aPd  0.288  to  O.329  m  for  sites  1, 

2,  and  3,  respectively. 


PART  IV:  CONCLUSIONS  AND  RECOMMENDATIONS 


Conclusions 


71.  Based  on  the  test  results  reported  herein,  the  following  con¬ 
clusions  are  considered  warranted: 

a.  Swept -frequency  radar  measurements  over  the  frequency 
range  of  0.5  to  2.0  GHz  can  be  used  to  estimate  power 
reflectance  from  the  surface  material  of  highway  pave¬ 
ment  structures  and  to  determine  the  amplitude  of  the 
subsurface  contribution.  These  measurements  make  pos¬ 
sible  estimates  of  the  electrical  properties  of  each  of 
the  materials. 

b.  Layered  materials  produce  interference  patterns  in  the 
power  reflectance  curves  over  the  frequency  range  of 
0.25  to  2.0  GHz.  Interference  patterns  can  be  used  to 
calculate  both  the  optical  thickness  of  each  layer  and, 
with  an  adjustment  made  for  wave  velocity  in  the  mate¬ 
rial,  the  physical  thickness  of  each  layer. 

£.  Insufficient  data  were  obtained  in  this  program  to  show 
whether  the  frequencies  of  2.0  to  8.0  GHz  are  as  effec¬ 
tive  as  those  of  0.25  to  2.0  GHz  in  determining  layered 
material  properties. 

Recommendations 


72.  It  is  recommended  that: 

a.  Field  tests  be  conducted  on  portland-cement  concrete 
highways  to  determine  the  capability  of  the  swept- 
frequency  radar  system  to  measure  depths  to  reinforcing 
steel  and  layer  thicknesses  in  the  pavement  structure. 

b.  Additional  tests  be  made  at  frequencies  of  2.0  to  8.0  GHz 
to  determine  the  effectiveness  in  measuring  surface  and 
subsurface  material  properties  at  these  ranges.  This 
would  require  the  procurement  of  more  efficient  antennas 
to  propagate  the  energy  at  these  frequencies  and  the  mod¬ 
ification  of  the  signal  recordi rg  equipment  to  increase 
the  accuracy  of  the  data. 

c.  The  system  be  modified  (both  the  system  size  and  method 
of  recording)  so  that  continuous  records  can  be  made 

of  layer  thicknesses  as  a  vehicle  containing  a  swept- 
frequency  radar  system  travels  along  the  highway. 
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I.r\ 


Table  1 


Highway  Data 


Site  1 

Site  2 

Site  3 

(Sta  724) 

(Sta  654) 

(Sta  604) 

Measured  asphaltic  concrete  thickness 

in  traffic  lane  X^  ,  m 

Measurements  at  center  of  test  site: 

0.226 

0.210 

0.214 

Asphaltic  concrete  thickness  at 

edge  of  shoulder  Xp  ,  m 
Clay-gravel  base  thickness  at 

0.146 

0.124 

0.102 

edge  of  shoulder  X3  ,  m 

Asphaltic  concrete  thickness  at 

0.286 

0.359 

0.483 

edge  of  emergency  lane  Xg  ,  m 
Clay-gravel  base  thickness  at 

0.121 

0.098 

0.083 

edge  of  emergency  lane  X^  ,  m 

Measurements  I.83  m  east  of  test 
site  center: 

0.311 

0.394 

0.552 

Asphaltic  concrete  thickness  at 

edge  of  shoulder  Xp  ,  m 
Clay-gravel  base  thickness  at 

0.121 

0.121 

0.114 

edge  of  shoulder  X3  ,  m 

Measurements  I.83  m  west  of  test 
site  center: 

0.305 

0.375 

0.470 

Asphaltic  concrete  thickness  at 

edge  of  shoulder  Xp  ,  m 
Clay-gravel  base  thickness  at 

0.140 

0.130 

0.095 

edge  of  shoulder  X3  ,  m 

0.305 

0.371 

0.489 

Maximum  thickness  of  clay-gravel 

base  in  traffic  lane  X5  ,  m 

0.186 

0.262 

O.362 

Minimum  thickness  of  clay-gravel 

base  in  traffic  lane  Xg  ,  m 

Moisture  content  measured  during 
highway  construction,  $: 

0.133 

0.209 

O.309 

Clay-gravel  base* 

4.4 

4.5 

4.8 

Silt  foundation 

Dry  density  measured  during 
highway  construction,  g/cnP: 

22.0 

18.0 

16.6 

Clay-gravel  base 

2.04 

1.94 

2.19 

Silt  foundation 

1.64 

1.58 

1.64 

*  Average  moisture  contents  of  clay-gravel  base  material,  measured 
shortly  after  radar  measurements  were  made,  were  7-6%,  7.2$,  and  6.8$ 
for  sites  1,  2,  and  3>  respectively.  The  moisture  content  samples  were 
taken  at  the  edge  of  the  asphaltic  concrete  near  the  bottom  of  the  clay 
gravel  base. 
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Table  2 

Results  From  Highway  Data  Obtained  with  Syre^t-Fre^uenc; 

(0.5  to  2.0  GHk) 
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FREQUENCY,  OH* 


PLATE  2 


0.5 


FREQUENCY,  CM* 


0.3 


^S^53fS£eK 


0.4 
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PLATE  7 


FREQUENCY 


PLATE  9 


CORRELATION  COEFFICIENT 


■  1  I  .  AWW'\A/>A/V,/AW>',‘,mMv/Av 

0  1.0  2.0  3.0  4.0  8.0  0.0  7. 

TEST  I 


TEST  2 


0.010  f-  1  f. 


1.0  2.0  3.0  4.0  3.0  3.0  7. 

OPTICAL  DEPTH,  IP 

TEST  3 


OPTICAL  THICKNESS 
HIGHWAY  TESTS 
SITE  I  (STA  724) 
TESTS  I,  2,  AND  3 


CORRELATION  COEFFICIENT 


TEST  4 
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)  1.0  2.0  3.0  4.0  3.0  0.0  7.0 

TEST  5 


0.003  i-i  •'  k 


^  il  ^  j  i  Kv \v.'VvV  X  •v'A/Vv'\'V>  v  V\ 

1.0  2.0  9.0  4.0  3.0  0.0  7.0 

OPTICAL  DEPTH,  M 

TEST  6 


OPTICAL  THICKNESS 
HIGHWAY  TESTS 
SITE  2  CSTA  654) 
TESTS  4,  5,  AND  6 


APPENDIX  A:  FOURIER  TRANSFORM  ANALYSIS 


2.  As  shown  in  fig.  Al,  the  points  defined  in  the  frequency  do¬ 
main  will  be  l/T  apart. 

3.  If  fix )  is  redefined  so  that  the  amplitude  from  to 

2T^  is  zero  and  the  new  period  T^  is  twice  as  long  as  ,  the 
spacings  between  the  points  in  the  frequency  domain  will  be  reduced  by 
half,  and  the  amplitude  of  C n  for  the  points  at  l/T^  >  2/T^  ,  3/T^  , 

...  frequency  will  be  reduced  by  half  (fig.  A2). 

4.  This  reduction  can  be  seen  in  equations  A3  and  A4  for  the 
Fourier  coefficients.  For  example,  from  equation  A3  for  a  period  : 

*  Refers  to  corresponding  listing  in  Literature  Cited  section  following 
main  text. 
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Fig.  A2.  Periodic  wave  and  frequency  function  after 
adding  zero-amplitude  points 
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5.  Since  the  function  /  (x)  is  defined  as  zero  between 

2T^  ,  the  second  integral  for  A^  in  equation  A?  Is  zero.  The 

cients  A  _  will  be  half  of  A  when  n_/2  is  enual  to  n, 

n2  nl  2  1 

of  the  2/T  factor  ir.  front  of  the  integral. 
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6.  This  same  process  can  be  generalized  for  the  case  of  a  func¬ 
tion  with  an  arbitrary  fraction  of  its  period  at  zero  amplitude.  The 
coefficients  for  the  nonzero  function  can  be  written  as 


Ani  “  T 


A 


'n2 


(A8) 


for 


T1 

ni  =  n2  t: 


Afll  for  other  than  integral 
n^  have  some  error  because  the  function  /^(x)  can  go  to 


7.  The  Fourier  coefficients  for 
values  of 

zero  in  the  middle  of  a  sine  or  cosine  cycle.  This  effect  can  be  re¬ 
duced  by  subtracting  the  mean  of  /- (x)  from  each  data  point  in  /  (x) 

A.  rj  X 

before  calculating  the  Fourier  transform.  This  has  no  effect  on  the 

A^^  values  for  integral  values  of  n^  .  Also  note  that  the  values  for 

A  for  fractional  values  of  n  less  than  1  have  no  meaning  because 
nl 

the  original  data  record  /^(x)  would  not  be  as  long  as  the  sine  or 
cosine  period  of  the  Fourier  transform. 
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